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ABSTRACT 


The  applicaton  of  simulation  techniques  provides  a  convenient  and  useful 
means  for  studying  and  developing  signature  analysis  schemes  to  describe 
orbiting  body  configuration,  orientation  and  motion. 

This  paper  discusses  a  simulation  of  the  overall  body  motion  characterized 
by  the  orbital  and  local  motion  about  the  center  of  mass.  This  motion  is  rela  ¬ 
ted  to  the  observations  made  by  ground  stations  and  provides  orientation  of  the 
total  motion  in  the  antenna  axes  systems  and  orientation  of  the  lines  of  sight 
relative  to  the  ground. 

Since  only  the  basic  characteristic  motions  are  required  and  exacting 
simulation  to  some  existing  satellite  is  not  necessary,  certain  simplifications 
are  introduced.  These  include  omission  of  perturbations  and  observational 
errors  and  the  elimination  of  a  true  time  requirement  and  a  time  based 
progression  reference. 

The  emphasis  here  is  on  a  description  of  the  equations,  relations,  and 
transformations  of  the  simulation  model.  The  computer  program  based  on 
this  simulation  model  is  described  in  a  separate  report. 
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1.0  INTRODUCTION 


The  overall  motion  of  a  satellite  is  characterized  by  its 
orbital  motion  and  local  mction  about  its  center  of  mass.  It  is  of 
interest  to  discern  both  these  motions  from  ground  viewing  stations. 

Thus  it  is  desired  to  relate  the  motion  to  observations  at  a  complex 
of  earth  stations  each  of  which  totally  positions  the  body  in  space 
relative  to  its  line  of  sight 

For  a  given  system  and  operation  mode  correlation  of  the  observa¬ 
tions  to  the  body  motion  can  provide  a  method  for  developing  signature 
analysis  techniques  for  describing  body  configuration,  orientation, 
and  motion.  Simulation  of  the  body  motion  and  measuring  techniques 
provides  a  convenient  approach  for  carrying  out  such  a  study. 

This  paper  will  be  restricted  to  a  discussion  of  the  simulation 
of  the  body  motion  only.  It  will  relate  this  motion  as  a  final  output 
to  a  receiving  antenna  axic  system  along  the  line  of  sight.  What 
is  discussed  for  one  site  carries  over  to  any  other  site. 

The  equations  relations,  and  transformations  required  in  the 
consideration  of  overall  motion  have  been  considered  in  earlier  internal 
reports  by  the  author.  The  latter  reference  includes  a  simulation 
of  the  resulting  scattering  matrix  from  stored  scattering  data  and  sets 
up  radar  outputs  for  variou-  polarization  modes.  This  radar  simulation 
also  includes  the  effect  of  Faraday  rotation  and  noise. 

In  a  computer  program  many  variations  are  possible.  These  include 
differences  in  describing  the  station >  orbit,  and  their  combination 
for  observation.  The  present  paper  presents  a  somewhat  simplified  ver¬ 
sion,  as  compared  to  earlier  reports  by  the  author,  of  the  motion  simulation. 

Since  only  the  ba^ic  motions  are  required  no  perturbations 
or  observational  errors  are  included.  Again,  since  it  is  not  the 
intention  here  to  simulate  motion  of  any  existing  satellite  as  close 
as  possible,  accounting  for  true  time,  taking  time  as  a  progression 
reference,  and  solving  the  transcendental  Kepler  equation  (see  section 


*  A  description  of  the  program  already  completed,  will  be  given  in  a 
subsequent  report. 
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on  orbital  motion,;  also  become  unnecessary. 

Thus,  using  a  revolving  spherical  earth  and  earth  surface  fixed 
site  locations  (earlier  internal  reports  of  the  author  comment  on 
ellipsoidal  earth  and  nonzero  altitude  stations),  the  orientation  of 
body  fixed  coordinates  with  respect  to  antenna  coordinates  are  deter¬ 
mined  with  time  Tbe^e  determinations  are  specified  in  terms  of 
three  orientation  angles.  in  addition,  the  orientation  of  the  station 
line  of  sight  relative  to  the  ground  arid  the  station  to  body  range 
are  provided  Iri  the  r^dar  case,  these  angles  are  used  to  look  up 
stored  scattering  data  obtained  either  from  measurement  or  calculation, 
Both  stabilized  and  unseabi 1 i zed  motion  of  the  body  are  considered. 
It  is  in  the  latter  no  torque,  case  that  the  simulation  of  local  motion 
enterb.  The  stabilized  modes  of  major  importance,  spin  stabilized, 
horizon  stabiLized^  and  earth  centered  stabilized,  are  accounted  for 
by  establishing  the  proper  orientation. 


2.0  THE  SIMULATION  MODEL 

The  basic  technique  of  trie  Simula:  ion  involves  a  sequence  of  six 
coordinate  transformations.  Jhe  coordinate  systems  dealt  with  are 
designated  in  order  as  follows. 


1 , 

geometrical  bodv  a*e.b, 

4 

describes  body  geometry  and 
svmme  tr le  s  t 

2, 

inertial  axes 

!bj 

p ,  describes  conveniently  body  motion 
about  center  of  mass 

J, 

angular  momentum, 

used  as  a  reference  with  respect 
tc  body  local  motion 

4 

orbital 

{'j 

^  used  as  a  reference  with  respect 
to  earth  motion 

5. 

geocentric , 

(gJ 

1*.  used  as  a  reference  with  respect 
to  earth  motion 

6 , 

topocentr ic , 

jsf,  moves  with  earth  and  locates  the 
'  J  site (s) 

7. 

radar , 

1  ! 
lrJ 

moves  with  radar  and  accounts  for 

antenna  orientation 


t  In  the  radar  case,  this  provides  a  reference  to  which  scattering  data 

is  related  in  terms  of  a  radar  measurement  coordinate  system  H-  Details 

are  given  in  an  internal  report  by  the  author  as  well  as  in  reference  1. 
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With  perturbations  on  orbital  and  local  motion  due  to,  for  example, 
earth  oblateness,  drag,  solar  system  influences,  solar  radiation,  and  gravity- 
magnetic  torques  omitted,  the  jgj*,  and  systems  are  inertial,  t 

Also  a  Keplerian  orbit  in  a  true  square  law  force  field  is  used.  In 
calculating  the  earth  rotation  rate  ^ ,  account  must  be  taken  of  the 
earth  revolution  which  effectively  increases  the  rotation  of  the  4 
frame  with  respect  to  the  */gj-  frame,  for  existing  satellite  simulation. 

There  are  six  linear  transformations  relating  each  coordinate 
system  with  the  one  which  follows  it  in  the  list  above.  These  trans¬ 
formations  are  thus  designated: 

Tb,b’V?’  T?,§’  Tr,g’  Tg,s’  TS,r 

The  relation  of  the  |b  J*  system  to  the  ^b^  system  can  be  specified 
calculated  from  body  mass  data  from  diagonal izat ion  of  the  inertia 


or 


tensor  . 


%  T 


and  I  are  time  varying  where  T  y  carries  the 
s  ,  r  h  ,  _ 

local  motion,  T  carries  earth  rotation  and  T  reflects  orbital 
g,s  s,r  f  ^ 

motion  which  is  specified  by  the  radius  vector  RP  in  the  \ r f  system. 

r  v  ^ 

Local  motion  is  specified  by  movement  of  the  ibf  inertial  axis  system 

V  L  J 

with  respect  to  the  i'f  system.  The  overall  orientation  of  the  body 

r  ^  f  ? 

axes  *jbfv  to  the  radar  system  vr  r  is  given  by  the  product  of  linear 
transformations  from  which  three  orientation  angles  are  determined. 


T,  ,  =  T  T  T-  IV  T  y  Tu 

b'r  s,r  g,s  §,g  5,5  b ,5  Mb 


A  brief  discussion  related  to  the  different  transformations  is 
appropriate  prior  to  presenting  the  overall  simulation  logic. 

Each  of  the  transformations  can  be  given  with  respect  to  a 

coordinate  system  as  a  3  x  3  matrix,  T  relating  the  a  system  to 

CX ,  p 

the  P  system.  More  precisely,  T  converts  the  coordinates  of  a 

f  \  a  ’ :  ^ 

vector  in  the  jaj  to  its  coordinates  in  the  -jpj  system.  on 

the  other  hand,  relates  the  basis  vectors  of  the  ■jgj'  system  from  the 
basis  vectors  of  the  ■faj*  system. 


Additional  detail  is  given  in  an  internal  report  by  the  author. 
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in  terms  of  the 


t 

In  general 


a, ^ 


■  ‘W 


It  is  sometimes  easier  to  consider  the  T 

a, 3 

T  g  transformations  as  successive  axes  movements. 

Certain  of  the  T  ^  are  considered  composed  of  factor  trans- 
formations  each  of  which  are  pure  rotations  (in  some  plane)  Tw , 
that  (T^)  '  =  (T^)”1  . 


so 


Then 


t 


Prime  is  transpose  and 


)  is  inverse 


4 


/  CXi 

Os 

\ 

Tb>  - 

02 

03 

\ 

Ya 

y3 

/ 

=  (Ctr1  CCg  Og')1 


<k  > 

0i, 


Yi  > 


OLg  ,  CC3  are  the  direction  cosines  of  onto 
@2 ,  03  are  the  direction  cosines  of  onto 
Y2 ,  y3  are  the  direction  cosines  of  onto 


The  Ix  ,  I2  ,  I3  signify  the  moments  of  inertia  along  the  correspond¬ 
ing  principal  moments  of  inertia  axes  ,  and  respectively. 

The  jb'|  system  has  its  coordinate  axes  along  directions  determined  by 
body  geometry  and  symmetries.  In  general  the  ^bj*  and  systems 

will  not  coincide.  The  specification  of  T*  ,  is  effected  in  a  number 

b ,  b 

of  ways. 

1.  By  inserting  as  input  data, 

2.  By  calculation  from  body  mass  and  configuration  data, 

3.  By  specification  as  identity  (T^,^  =  I). 

The  calculation,  if  performed,  is  done  so  only  in  the  torqueless, 
unstabilized  mode.  For  all  stabilized  modes  it  is  sufficient  to 
consider  no  difference  between  the  ^bj*  and  ^b'j-  systems  (T,  =  I) 

with  controlled  orientations  taken  with  respect  to  the  ^b'j  system. 


b,r 


CLi  is  a  column  vector  with  components  CLX  ,  ,  Yi  etc.  Indeed  CL{  ,  Otg  ,  ctj 

are  eigenvectors  of  the  moment  of  inertia  tensor. 
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is  the  direction  of  the  angular  momentum  vector  of  local  motion, 
1 ixed  in  inertial  space.  The  vectors  form  a  right  handed 

system  and  are  specified  with  respect  to  the  ^ f *  as  noted  later. 


Ihe  Luier  angles  9, 


measure  the  motion  about  the  center  of 


mass  of  the  body  inertia  axes  n  b  r  relative  to  the  fixed  in  space  system 
The  line  of  nodes  is  the  intersection  of  the  (X^,  Y^)  plane  and 
the  ,  §2 )  plane.  A  description  of  the  equation  describing  the  Eulei 
angles  as  functions  of  time  is  given  below  under  the  section  local 
motion.  We  have  v 


b,f 


V' 

^cos(-0) 


T  r  -  T<3~  (-<?) 
b  ,  .  b  ,  - 

T  (sM 

b,§ 

(-6) 

V 

sin(-C)  0  \ 

/  1 

0 

f 

y 

cos(-$)  0 

0 

cos(-9) 

o  J' 

\°  - 

•  s  in  (  ~  9  ) 

(-Y) 


sin(-Y) 

cos(-Y) 


cos (-9) 


0  0  \  /-sin  cos  f 

cos  9  -sin  9  1  I— cos  V  -sin  ’f 

sin  9  cos  9  /  \  0  0 


^  ^  (- f )  takes  Y  ) 
b ' 

fx; 

f2 )  *  J  k 

r  ’ (  -fi  )  takes  <  Y'  \ 

\  b 


into 


into 


x; 


y^h-h 


takes  { Y"\ 

,  b 


into 


t  In  general  the  angle  between  X^  and  ?i-  0+Y  as  b-0 . 
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For  spin  stabilized  mode,  i  =  is  fixed 

Y  =  0  (the  spin  rate  about  Z^)  is  fixed* 


in  the 


system  and 


Then  T 


b,? 


cos  Y  -sin  Y  0 

=|  sin  Y  cos  Y  0 

0  0  1 


(i.e. ,  0=0,  0=0,  Y  =  Yq  +  YAt) 

For  horizon  stabilized  mode,  Z,  9*  Z,  is  parallel  to  R-  and  thus  the 

b  b  £ 

system  is  bypassed  and 


{?} 


Tb,r  ■ 1 


V  V  V 

(i.e* ,  0  -  0  -  Y  -  0) 


For  earth  centered  stabilized  mode,  Z,  i  =  Z,  is  parallel  to  R»  and 

b  b  ^  ^ 

thus  similarly,  ■>  =  I  (i.e,,9=0  =  Y  =  O,' 


(not  in  plane  of  ) 
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Since  rotation  of  theCf^,  f2)plane  about  §3  is  arbitrary,  we  choose 

T-y  _  to  be  given  only  by  the  two  angles  ax  and  ctg  .  (In  fact,  <£(0), 

S  >  9  f 

which  measures  the  initial  rotation  of  the  *  b  f  system  about  f 3  ,  is 
specified  and  arbitrary  in  any  case.)  5:  is  in  the  direction  of 
perigee  from  the  <j5j*  origin  .earth  center)  >  £3  is  normal  to  the  orbit 
plane,  and  5a  f°rms  a  right  handed  coordinate  inertial  system. 


F 


» 


Gfi  > 


~~  ( 1  ) 

Ty  }  i~a. ) 


COs 

Os 

-sin 

°\ 

/  COS  Cli 

0 

-sin 

ain 

cos 

Oj, 

0 

1  0 

1 

0 

0 

0 

1/ 

\  sin  CT: 

0 

cos 

«l/ 

For  unstabilized  motion  &i  and  a2  are  specified  as  input. 

For  spin  stabilized  motion  C.±  and  Og  are  specified  as  input, 

For  horizon  staDilized  motion  a*  *  TT/2  is  specified  and  Ug-  02 

(places  Z  pin  parallel  with  R-) 
b  5 

For  earth  centered  stabilized  0^  &  tt/  2  is  specified  and 

(^places  Z.  ,  in  parallel  with  R-) 

The  time  variation  of  the  radius  vector  of  the  orbiting  body,  R_  , 
and  the  de termination  of  the  angles  0i  and  £aare  given  below  in  the 
section  on  orbital  motion,, 


Z 

*  g 

I 


a 


The  geocentric  coordinate  system  /gj’  is  chosen  in  an  arbitrary  way 

since  initial  time  is  of  no  consequence  for  the  purpose  of  this 

simulation.  In  other  words,  the  X  direction  is  taken  as  the  direction 

g 

in  space  made  by  the  earth  radius  passing  through  the  equator  at  0° 

longitude  at  an  arbitrary  initial  time,  t  ,  which  is  taken  as  zero. 

The  Zg  is  the  direction  from  earth  center  through  the  north  pole  at  t  . 

Y  completes  a  right  handed  inertial  system.  With  no  perturbations 
g 

Cl  are  time  invariant  where: 

\d  =  inclination 

A  =  argument  of  perigee  t 

i-  -  right  ascension  of  ascending  node 


T?  =  Tp  = 
?  g  §  >  g 


COS 

-sin  L 

0\ 

/ 1 

0 

0 

\ 

f  cos 

A 

-sin 

A 

°\ 

sin 

cos  C 

0 

1° 

cos 

-sin 

sin 

A 

cos 

A 

0 

0 

0 

1/ 

\0 

sin 

M- 

cos 

0 

0 

i/ 

g>s 


Perigee  or  apogee  are  convenient  choices  of  orbit  epochs  in  an  orbit 
period . 


and 
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The  2^  direction  is  normal  to  earth  surface  pointing  away  from  the  earth 

at  the  site  while  X  and  Y  determine  the  tangent  plane  at  the  site 

s  s 

with  X  along  a  meridian  to  the  north  and  pointing  west  where 

•  =  longitude  of  site  \  o  *■  X  2T  ) 

r\  =  measures  latitude  of  site  (-tt/2  s.  x\  £  tt/  11) 


j  cos  p 

0 

sin  3  ^ 

l 

/  COS 

1 

a 

s  in 

a  0  \ 

t  =  i 

0 

1 

0 

j  I 

-s  in 

CL 

cos 

a  0 

s-s  i 

^-sin  3 

0 

cos  p  j 

1  1 

l  0 

0 

W 

where  a 

+  X 

+  TT  ; 

At 

=  t 

-  co 

(with  tQ  arbitrary  we  usually 


take 


=  7]  -  TT/2 


T 

s,r 


t With  simulation  of  actual  orbits,  fixing  ^gj  in  space  at  tQ  consistent 
with  astronomical  data  would  mean,  for  example,  taking  tQ  as  vernal 
equinox  time. 
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R  is  the  radius  (line  of  sight)  vector  to  the  orbiting  body  expre  w 

^  f  )  /\  A 

in  the  s  j  system,  and  the  look  angles  9  and  0  will  be  cur^i  k  . 

in  the  section  on  orbital  motion. 

The  Y  of  the  radar  system  lies  in  the  (Xg  ,  Ys)  plane,  Lhat  i. 

A 

remains  horizontal.  It  then  is  the  axis  about  which  elevation,  9> 
changed,  $  measures  azimuth  change. 


{  cos 

a 

6 

0 

sin 

/  cos 

A 

0 

A 

sin  0 

0  \ 

0 

1 

0 

-sin 

$ 

A 

COS  0 

0 

(-sin 

A 

6 

0 

cos 

S) 

1  o 

0 

1  / 

T  ,  =  T 
b,r  s  , r 


^  In  the  radar  case,  the  line  of  nodes  coincides  with  the  ,  of  the 
radar  measuring  system  jr'J*  which  is  a  rotation  of  jr J  about  Zr  of 
-03 ,  Yr  and  Yr ‘  are  horizontal  axes  in  their  respective  systems 
9  and  0  are  drawn  to  be  consistent  with  an  adopted  use  in  the  |r 1 \ 
system  for  scattering  measurements  (see  reference  1  ) . 
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al  1 

ala 

al  3 

a2  1 

a22 

a23 

a:H 

a3  2 

a3  3 

I 

f  cos  P 3 

sin 

/  1 

0 

0 

N  jCOS 

3a 

sin  32 

-( 

-sin  p3 

cos  p3 

0 

11° 

cos  ^  x 

sin 

3i  lj-sin 

32 

cos 

° 

[  0 

0 

1  / 

\  0 

\ 

-sin  Pi 

cos 

0 

0 

This  gives 


a:  l 

-  cos 

-*> 

p3 

cos 

3  2 

-  sin 

33  cos 

3i 

s  in 

02 

d2 1 

--sin 

3 

3 

cos 

32  ■ 

-  cos 

33  cos 

0i 

sin 

02 

a3  1 

=  sin 

Si 

s  in 

32 

al  2 

=  cos 

3 

3 

s  in 

0 

2 

4  s  in 

p  cos 

3 

01 

cos 

32 

a2  2 

=-  sin 

33 

s  in 

4-cos 

3-3  cos 

0i 

cos 

3  2 

a3  2 

=  -  s  in 

0i 

cos 

al3 

=  sin 

p  3 

s  in 

3i 

a23 

=  cos 

3'3 

sin 

0i 

a3  3 

=  cos 

^  1 

Then 

-l 

(1 

-43) 

12 

1 

1 

Pi  = 

can 

a3  3 

—  LUS 

a3  3 

q 

-2 


-1 


tan 


c 


a3  1  "N 
-%S  ' 


rsZ 

0  -  ?!  ,-0  =  (3  n/2  -  P  ),  0  =  P  2 

6*,  0,  and  p3  are  the  final  antenna  body  orientation  angles. 


The  overall  simulation  logic  can  now  Lc  conveniently  displayed  schemati¬ 
cally  by  the  following  diagram*  Of  course  either  stabilized  motion  or 
local  motion  is  appropriate  For  convenience  both  are  shown.  The 
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symbol  (x)  designates  alternate  paths.  The  logic  flow  chart  for  the 
actual  program  appears  in  a  subsequent  report. 


3  0  ORBITAL  MOTION 


With  an  inverse  square  law  radial  force  field  from  earth  center, 
the  equation  for  motion  in  the  system  becomes 


R§ 

R? +  d  jC  = 


=  0 


with  d  =  constant,  accounting  for  gravity  and  the  mass  of  the  two 
body  system.  In  the  metric  system,  for  example 

d  =  3.9858  x  105  km3/sec? 

This  results  in  planar  motion,  that  is  R  x  R  =  constant  so  that 


R  =  a  R  +  b  R  :  with  a  and  b  scalar  functions  of  time  and 
o  o 


R  and  R  evaluations  at  some  time  t1. 
o  o  o 

Choosing  t^  to  be  the  time  of  perigee  and  the  coordinate  system 

then  places  R^  along  5i  >  a^-onS  ?s  an<^  gives 

2  \ 

Rc(t)  =  a(cos  E-e)  51  +  a(l-e  )  sin  E 


and 


R.(t)  =  -a  E  sin  E  +  a(l-e2;  E  cos  E 


where  a  is  the  semimajor  axis  of  the  elliptical  orbit  and  e  is  the 
ellipticity  of  the  orbit.  E(t),  the  eccentric  anomoly 
is  related  to  t  in  Kepler's  equation 


E-e  sin  E  -  (d/aJ)  At  ;  At  =  t  - 


(E  =  (d/a3)'  /  (1-e  cos  E)  ) 


14 


where  is  time  of  perigee  crossing,,  If  initial  time  tQ  is  not  t* 
the  Kepler  equation  is  slightly  more  complex  but  straightforward.  We 
take  tQ  =  t^.  Also  equal  increments  in  At  produce  unequal  increments 
in  E  and  vice  versa.  Solving  for  E  from  t  is  difficult  whereas  solving 
for  t  from  E  is  direct.  The  programmed  version  discussed  in  a 
subsequent  working  paper  uses  equal  increments  in  E  as  reference. 

For  simulations  of  actual  orbits  where  motion  is  referenced  to  real 
times,  solution  of  the  Kepler  equation  for  E  can  be  accomplished  using 
an  iterative  scheme  such  as  Newton-Raphson.+ 

From  the  expressions  for  Fp.  and  Rg.  we  have 


Si 

per igee 


■ 

Let  f (E)  =  (E-e  sin  E)  -A  ;  A(t)  =  fd/a3 ) 2 ( t- tQ) 


f (E)  =  df/dE  =  1-e  cos  E 

With  En  the  nc^  estimate  of  true  E,  let  hn  be  such  that 

f  (E  ) 

f(En)  4-  hn  f(En)  =  0,  i.e.,  hn  =  -  -r-2- 

f(En) 

We  take  En  +  ^  *  hR  +  En  and  repeat  the  process  for  hn+^ 
|f(En)|  <  e,  a  preselected  small  number. 
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until 


Nov  R  s'c.>  =  F  T_  R_  ; 

-  %  b  t  ,  V  ^ 

where  as  noted  before ,  p  is  earth 
radius  taken  away  from  earth  center. 


n  =  Rs  •  Xs  etc. 


cue  discussion  of  T  in  the 

s  ,  r 

pt  lUUt)  Sect  j  Ji:  „ 

A 

Of  various  acquisition  test-  a  useful  and  simple  scheme  compares  9 

A 

against  some  0X  such  that 

A  A 

if  'i)  0  <  0X ,  target  visible 

A  A 

(ii)  0  ^  0X  ,  target  not  visible 

If  (ii)  holds , proceed  along  orbit  until  (i)  holds. 

If  (i)  holds  .continue  with  main  calculations. 


-  x 


-2  .  +  r-3ZS 


?  \ 


the- 


A  l 

9  ^  t  a- 1 


r 


(  --  ) 


1 1 ie  c e oaiw  cry  or  9  and  0  v a s h o wo  in 


4.0  LOCAL  MOTION 

The  movement  of  any  vector  fixed  in  the  body  and  so  moving  with 
the  ^bj*  system  with  respect  to  the  system  can  be  described  by  the 

three  Euler  angles  0(t),  ¥(t),  and  0(t)  defined  in  the  discussion  of 


Indeed,  these  angles  satisfy  the  following  Euler  differential 


equations,  for  zero  applied  torque. 


We  call  the  plunged  condition  of  the  radar  with  0  +  TT  rather  than  0 
(also  possible ) .  Then  for  Xr  above  the  (Xg ,  Y Q)  plane  in  the  normal 
condition,  Xr  and  Rs  have  projections  in  the  plane  in  opposite  direc¬ 
tions  whereas  in  the  plunged  condition  Xf  is  below  the  (Xg ,  Yg)  plane 
and  the  two  projections  lie  in  the  same  direction. 
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9 


Ui-Ia) 

Ii  la 


=  M  sin  9  cos  Y  sin  Y  =  nutation  rate  (ioi. 

A9  variations) 


S; 


M  2  M  2 

0=7  cos  Y  +  —  sin  Y  =  precession  rate 
Li  J-a 

i  /  M  M  2wM.2w'N  ~ 

Y  =  v  7  -  7  cos  t  -  7  sm  Y  )  cos  9  =  spin  rate 

^  is  il  S 

where  for  generality^  we  take  Ii  ^  I2  ^  I3  as  the  three  moments  of 
inertia  with  respect  to  the  principal  axes  of  inertia  X^  ,  ,  Z  and 

M  is  the  magnitude  of  the  angular  momentum,  M  along  . 

Denoting  the  corresponding  angular  velocities  as  Qx  (t) ,  ( L )  anc 

fi3(t)  with  directions  along  the  X^ ,  Y^  ,  axes  we  have 


9(0)  =  cos 


k  0,  (Q)  J  ^  y 


M 


(0) 


and  M  =  T  I . 


i  =1 


i  i 


H3  =  l\  Ci\ 

i  =  l 


tan 


(0)  | 
Ii  Cii  (0)  J 


0(0)  is  arbitrary  and  is  specified.  Note  initial  values  of  9  ani 
reflect  initial  values  of  9,  0,  and  Y. 


The  solutions  to  the  Euler  equations  (see  reference  2O  are  obtained 

in  series  which  allow  for  approximations  to  any  desired  accuracy . 

ft 

These  solutions  are. *  11 


In  all  v-^ses  each  moment  of  inertia  along  an  axis  is  less  i/ha 
equal  to  the  sum  of  the  moments  of  inertia  along  the  other  rwo 
orthogonal  axes.  For  definetness  the  inequalities  are  used. 

1 1 

1  In  order  that  the  expressions  for  9(t),  Y(t),  and  0(t)  be  consist  nt 
with  initial  conditions,  the  t  used  in  these  expressions  should  be 
shifted  with  respect  to  true  (orbit)  time. 
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0(0 


-i  f(sinh  y-q  sinh  3v  +• 

co«  - "“a - - 

Hcosh  V+q  cosh  3v  +• 


»)(l+2q  cos  2n.t+2q4cos  4a t  +••») 

•)(l-2q  cos  2}Jt+2q4cos  4|Jt  -•••)  J 


from  cos  0 


£  (ia/2K)  ^(Xt/2K) 

i^10(ia/2K)  ^oi  (WW 


(^xx  O  are  the  four  types  of  theta  functions  and  K  is  the  complete 
elliptic  integral  of  the  first  kind.) 

,  _  “l  f(l+2q  cosh  2y-+-2q4cosh  4y  H-***)(sin  ^t-q2sin  3|Jt+«»-)l 

1(1- 2q  cosh  2y-t-2q  cosh  4v -  •  •  • )  (cos  jJt+q  cos  3jat+*-«)  J 


from  tan  ¥  - 


2y+2q  cosh  4y - )  (cos  jJt+q  cos  3|at+*-.) 

^00(.ia/2K)  ^  (Xt/2K) 

JSoi  <  la/2K)  $io  (W2K) 


•2i0(t)=  Fit)  exp  f 2i  [-  +  4a  (  q  sinh  ,2v-2q4sinh  Uy  +•  • ’ 

•  p  l  j x  (  1  - 2q  cosh  2y+2q  cosh  Ay  +• 

X  oi  (W2K) 


T)]t} 


=  F(t)  exp  |2i  [  “ 


2iK  $  D1  (ia/2K) 


where 


_i-  K  *  /K 

e  ;  modulus  of  K  is  k. 


q=k°  /  1 6  +  k  /  3  2  +  -  -  •  ;  q  2  = 

(k'J*  -  (  c  «  i  1,3 

la'!-?  (  I,  C-n  )  1  1 


i=l 


=  modulus  of  K 1  =  1  -  k? 

a  ( I2  ~  I3  )  '.Cl 3.  -Ms  ) 

^  2K  ’  li  I3 13 

y  =  rra/K ,  (~)^  =  1  +  2q  +  2q4  +  2q9  +  ••• 

f  Ig  (!].  c-m2  )  "j  =  £00  (  ia/*0 
h;  (IjjC-M2  )J  C)01(ia/2K) 


F(t) 


B 


-V 
x  < 


Xt-  ia 
2K 

Xt+ia 

2K 


) 

) 


B  is  a  constant  specified  by  0(0). 
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Together  with  0(0),  q,  ja,  and  y  represent  a  possibly  more  convenient 
input  set  than  1^,  i=l,30  Indeed  for  q  «  1,  a  meaningful  associu 

tion  gives  q  as  a  measure  of  nutation  amplitude,  A9,  jj  as  a  measure  of 
spin  rate,  Y,  and  y  as  a  measure  of  precession  angle,  0.  Writing 


g2i0(t)_  e2i(0x  +  08 ) 


with  e 


2i0i 


F(t), 


0S 


_x_  ^<la/2K> 

2iK  201(la/2K)  ' 


<3'0  -  ^  ) 

then  02(t)  ==  constant,  measures  average  rate  of  0  variation  (precession 
rate)  and  0x(t)  accounts  for  (periodic)  variations  about  the  average 

0  (“0p)« 


The  frequency  of  the  0  and  Y  periodicities  = 
The  frequency  of  the  0X  component  of  0  =  ^, 


The  frequency  of  the  02  component  of  0  B  TT0a  . 


Various  interesting  special  cases  occur.  For  example  constant 
precession  (0X  =  0)  is  obtained  about  §3  with  spin  about  when  1!=  I2  . 
(See  reference  1.)  For  unstabilized  bodies  the  particular  motions  of 
interest  are  normally  tumbling,  precession  and  nutation.  In  the  long  (ti.k 
term,  the  dominant  motion  is  usually  tumbling.  All  of  these  motions 
represent  simplifications  of  the  general  motion. 
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